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We present results from an analysis of B° p~ decays using (383.6 ± 4.2)xlO®BS pairs 

collected by the BABAR detector at the PEP-II asymmetric-energy B Factory at SLAC. The mear 
surements of the — » p~ branching fraction, longitudinal polarization fraction /i, and the 
CP-violating parameters Siong and Ciong are: 

B{B°^p+p-) = (25.5±2.1(stat)^^:^(syst))xl0"^ 

fL = 0.992 ±0.024(stat)+°;°i3(syst), 

Siong = -0.17±0.20(stat)+°;°^(syst), 

Ciong = 0.01 ± 0.15(stat) ± 0.06(syst). 

We determine the unitarity triangle angle a, using an isospin analysis of B — > pp decays. One 
of the two solutions, a = [73.1, 117.0]° at68%CL is compatible with standard model-based fits of 
existing data. Constraints on the unitarity triangle are also evaluated using an SU{3) symmetry 
based approach. 
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I. INTRODUCTION 

Charge conjugation-parity (CP) violation was first 
seen in the decay of neutral kaons [1]. ft was shown 
some forty years ago that CP violation is a necessary but 
insufficient condition required to generate a net baryon 
anti-baryon asymmetry in the universe [2] . The standard 
model (SM) of electroweak interactions describes CP vi- 
olation as a consequence of a complex phase in the three- 
generation Cabibbo-Kobayashi-Maskawa (CKM) quark- 
mixing matrix [3, 4]: 

VcKM = Kd Kb . (1) 

V Vtd Vts Ftb / 

VcKM describes the couplings of the u, c and t quarks to 
d, s and b quarks, which are mediated by the exchange of 
a W boson. In B-meson decays the CP violating param- 
eters of the SM arc most directly related to the angles 
and sides of the so-called Unitarity Triangle (UT), shown 
in Fig. 1. The angles a, (3, and 7 are defined as 

a = arg[-14dnb/KdKrb], (2) 
I3=^vg[-V,^V:JV,AV:^], (3) 

7 = arg[-KdKrb/KdV;b]- (4) 

Any non-trivial phase in Vlj is CP violating. CP violat- 
ing phases originating from the CKM matrix are related 
to weak interactions, and therefore referred to as weak 
phases. In the Wolfenstein convention [5], the angle 7 
is the phase of Kib, P is the phase of T4d, and a is the 
phase difference between V^h and Vu\ constrained to sat- 
isfy a = TT — /3 — 7 through the unitarity of Vckm- 




(0,0) (1,0) 

FIG. 1: The Unitarity Triangle in the p-^ plane. 

In Fig. 1 the UT is shown in the complex (p, if) plane, 
where the apex is given by the phase- convention inde- 
pendent definition [6]: 

- , — _ KudV' b 
VcdVcb 
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The quest to understand CP violation remains, de- 
spite its observation in the B meson system by both the 
BABj^R [7] and Belle experiments [8], since the SM does 
not, through the CKM phase, incorporate enough CP 
violation to explain the current matter- antimatter asym- 
metry [9, 10]. The CKM angle p is measured to a 1° pre- 
cision with h CCS decays [11, 12] and is consistent with 
current predictions [13-15]. A significant deviation from 
results of SM-based fits of existing data for any of the UT 
angles would be a clear indication of new physics [16, 17]. 

In principle, direct experimental measurements of a 
can be made from decays that proceed mainly through a 
h uud tree diagram such as B° p+p", p^TT^, tt+tt" 
and a^TT [18, 19]. Interference between the direct de- 
cay and decay after B'^B^ mixing in each of these decays 
results in a time-dependent decay-rate asymmetry that 
is sensitive to the angle a. Figure 2 shows the leading 
order tree and gluonic penguin loop contributions to the 
decay 5° P~ ■ The presence of penguin contri- 

butions complicates the extraction of a. from these de- 
cays. Using isospin relations [20], measurements of the 
B+ ^ p+p° [21, 22] and B° p°p° [23] branching frac- 
tions show that the penguin contribution in B ^ pp is 
smaller than the leading tree diagram. The use of SU (3) 
flavor symmetry to increase the precision on the weak 
phase constrained using B" P^P~ and i?+ K*^ p'^ 
decays has also been proposed [24]. Both of these ap- 
proaches are discussed in Section X. Section X and 
Ref. [25] describe a number of possible sources of the- 
oretical uncertainty. 




d 



FIG. 2: Tree (top) and gluonic penguin (bottom) diagrams 
contributing to the process p~ ■ The penguin con- 

tribution coming from the diagram with a top quark in the 
loop dominates as contributions from processes with u and c 
quarks are suppressed. 



In P" decays, a spin particle (the B° meson) 

decays into two spin 1 particles [p^ mesons), as shown in 
Fig. 3. Subsequently each p^ meson decays into a tt^tt" 
pair. As a result, the CP analysis of B^ P'^ P~ is com- 
plicated by the presence of one amplitude with longitudi- 
nal polarization and two amplitudes with transverse po- 
larization. The longitudinal amplitude is CP-even, while 
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the transverse amplitudes contain CP-even and CP-odd 
states. The decay is observed to be dominated by the 
longitudinal polarization [26, 27], with a fraction fi de- 
fined as the fraction of the helicity zero state in the decay. 
Integrating over the angle between the p decay planes 
the angular decay rate is 



are defined as 



Td cos Old cos 62 4 



Jl cos^ 9i cos^ O2 



+ -(1-/L)sin2 Or sin2 62 



(6) 



where the helicity angles 0i=i,2 are the angles between 
the TT^ momentum and the direction opposite to that of 
the B'^ in the p rest frame. 




FIG. 3: A schematic of the decay of a B meson via two p 
mesons to a four pion final state. The p meson final states 
are shown in their rest frames, and is the angle between the 
decay planes of the p mesons. 



We identify (tag) the initial flavor of the reconstructed 
B candidate, B^^^i using information from the other B 
meson, iJtag, in the event. The decay rate f+{f-) for a 
neutral B meson decaying into a CP eigenstate, where 
the Stag is a B° is given by: 



/±(Ai) = 



-lAtl/Tso 
4tso 



±r]f[Ssm{AmdAt) - Ccos(AmdAt)] 



(7) 



Here At is the proper time difference between the decay 
of the Srec and Stag mesons in an event, tbo is the mean 
B° lifetime, Ama is the mixing frequency [28] , and 

?7/ is the CP eigenvalue of the decay. This assumes that 
there is no difference between the and lifetimes, 
AT = 0. The parameters S and C describe B-mixing 
induced and direct CP violation, respectively. For the 
longitudinal polarization 77/ = 5* = S'long and C = 



*^long — 



2 2m Xcp 
l + |Acp|2 



1- |A, 



■CP\ 



ong 



l + |Acp|2' 



(8) 
(9) 



where Xcp = ^ ^ [29], q and p are complex constants that 

relate the B meson flavor cigcnstatcs to the mass eigen- 
statcs, and A/ A is the ratio of amplitudes of the decay 
of a or B" to the final state under study. CP viola- 
tion is probed by studying the time-dependent decay-rate 
asymmetry 



A{At) 



r{At) - Y{At) 



(10) 



r(Ai) r(Ai) ' 

where r(r) is the decay-rate for B^ {B^) tagged events. 
This asymmetry has the form 

A(At) = Ssm(AmdAt) - C cos(AmdAt). (11) 

The transverse polarization is an admixture of CP-even 
iVf = +1) and CP-odd {rjf = —1) parts, where each CP 
eigenstate has a decay-rate distribution of the form of 
Eq. (7). As the longitudinal polarization dominates this 
decay, we describe the CP admixture of the transverse 
polarization with common parameters S'tran and Ctran- 
We set S'tran = Ctran = 0, and vary these parameters 
when evaluating systematic uncertainties. 

In the absence of penguin contributions in B° — !■ p'^p~ , 
•S'long = sin 2a and Ciong = 0. The presence of penguin 
contributions with different weak phases to the tree level 
amplitude shifts the experimentally measurable param- 
eter Q!eff away from the value of a. In the presence of 
penguin contributions acU = n + Sa, where 



'^long = - C^ong sin 2aeflF, 



(12) 



and Ciong c^an be non-zero. 

For _B+ P^ P^ one measures a time-integrated charge 
asymmetry, which is analogous to Eq. (10) and is defined 
as 

Acp = (13) 

N + N ^ ' 

where {N) is the number of B {B) decays to this final 
state. Table I summarizes the most recent measurements 
of the complete set of experimental inputs needed to ex- 
tract a from the B pp decays. 

In this article, we present an update of the time- 
dependent analysis of B^ P'^P~ and measurement of 
the CKM angle a reported in Ref. [26] and branching 
fraction reported in Ref. [30]. Improvements to the anal- 
ysis reported here include an increase in data sample an- 
alyzed, a tighter constraint on the proper time difference 
between the two B meson decays in selected events, an 
improved algorithm to determine the flavor of B mesons, 
a modified multivariate analyzer for background suppres- 
sion, and an improved description of the signal and back- 
ground model. 
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TABLE I: Measurements of the B pp decays. Branching fractions {B) are in units of 10 ®. The column marked Ciong(^cp) 
denotes the time dependent charge asymmetry Ciong in the case of the — » p'^ p~ decay, and the time integrated charge 



asymmetry Acp in the case of B 



Mode 



Expt. (luminosity fb 



B ( xlO" 



Clong(^Cp) 



Slo 





P'^P~ 


BABAR. [30] (82) 


B° 


-» p+p- 


BABAR, [26] (210) 


B° 


P'^P' 


Belle [27] (253) 


B" 


-> p+p' 


Belle [31] (535) 


i3± 


- P±P" 


H4B.4R, [21] (210) 


B± 




BcUc [22] (78) 


B° 




BABAR [23] (349) 



0.978 ± 0.014 
0.9411°:°^^ ± 0^030 



0.021 
0.029 



30 ± 4 ± 5 
22.8 ± 3.8l2;6 

16.8 ±2.2 ±2.3 
31.7±7.ll87 
1.07 ±0.33 ±0.19 0.87 ±0.13 ±0.04 



0.905 ± 0.042lo:'o27 
0.95 ±0.11 ±0.02 



-0.03 ± 0.18 ± 0.09 

-0.16 ±0.21 ±0.08" 
-0.12 ±0.13 ±0.10 
0.00 ±0.22 ±0.03 



-0.33 ±0.241 



0.19 ±0.30 ±0.08 



"'Belle Collaboration uses the opposite sign convention for Cion 



II. THE DATASET AND BABAR DETECTOR 

The results presented in this paper are based on data 
collected with the BABMi detector at the PEP-II asym- 
metric e+e" storage ring [32] operating at the Stanford 
Linear Accelerator Center. At PEP-II, 9.0 GeV elec- 
trons and 3.1 GeV positrons are collided at a center- 
of-mass energy of 10.58 GeV which corresponds to the 
mass of the T{AS) resonance. The asymmetric energies 
result in a boost from the laboratory to the ccntcr-of- 
mass (CM) frame of (3^ « 0.56. The dataset analyzed 
has an integrated luminosity of 349 fb~^ corresponding 
to (383.6 ± 4.2) xlQf^BS pairs recorded at the r(45) res- 
onance (on-peak). An additional 27.2 fb~^ of data were 
recorded about 40 MeV below this energy (off-peak) for 
the study of continuum background, where light quarks 
are produced in the process e+e" ^ qq {q = u, d, s, c). 

The BaBar detector is described in detail else- 
where [33]. Surrounding the interaction point is a five 
double-sided layer silicon vertex tracker (SVT) which 
measures the impact parameters of charged particle 
tracks in both the plane transverse to, and along the 
beam direction. A 40-layer drift chamber (DCH) sur- 
rounds the SVT and provides measurements of the mo- 
menta for charged particles. Both the SVT and DCH 
are surrounded by a solenoid magnet, that provides a 
1.5 T magnetic field. Charged hadron identification is 
achieved through measurements of particle energy-loss in 
the tracking system and the Cherenkov angle obtained 
from a detector of internally reflected Cherenkov light. 
A CsI(Tl) electromagnetic calorimeter (EMC) provides 
photon detection, electron identification, and ir^ recon- 
struction. Finally, the instrumented flux return of the 
magnet allows discrimination of muons from pions. For 
the most recent 134.0 fb~^ of data, a portion of the re- 
sistive plate chambers constituting the muon system has 
been replaced by limited streamer tubes [34-36] . 

We use a right-handed coordinate system with the z 

axis along the electron beam direction and the y axis up- 
ward, with the origin at the nominal beam interaction 



point. Unless otherwise stated, kinematic quantities are 
calculated in the laboratory rest frame. The other refer- 
ence frame which we commonly use is the CM frame of 
the colliding electrons and positrons. 

We use Monte Carlo (MC) simulated events generated 
using the GEANT4 [37] based BaBar simulation. 



III. RECONSTRUCTION OF B CANDIDATES 

A. Photon and reconstruction 

Photons are reconstructed from localized energy de- 
posits in the EMC that are not associated with a charged 
track. We require photon candidates to have an energy 
greater than 50McV, and a lateral shower proflle [38] 
to bo consistent with the photon hypothesis. We recon- 
struct neutral pions from combinations of two distinct 
photon candidates where the invariant 77 mass m-y^ sat- 
isfies 0.10 < < 0.16 GeV/c^. 



B. reconstruction 

We combine reconstructed 7r° mesons with charged 
tracks that are consistent with the tt^ hypothesis to form 
p* (;andidates. The invariant mass m.n.±.^o of the recon- 
structed p"^ is required to lie between 0.5 and 1.0 GeV/c^, 
to select events in the vicinity of the p resonance. We 
require that the helicity angle of each p meson satisfies 
—0.90 < cos^i < 0.98. This selection criteria suppresses 
continuum and B backgrounds. 



C. B" reconstruction and event selection 

We combine two oppositely charged p candidates to 
form the B candidate Brec- All particles in the rest of 

the event (ROE) arc combined to form the other B meson 
in the event Btag. In addition to the two charged tracks 
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in the -Brec, we require that there is at least one track in 
the Stag- 

In order to suppress potential backgrounds from 

e~ — > e+e", events, we require the second- 

to-zeroth Fox- Wolfram moment i?2 [39] of the event to 
be less than 0.98. Continuum events are the domi- 
nant background which is reduced by requiring the ab- 
solute value of the cosine of the angle between the B^^c 
thrust axis (TB) and that of the ROE (TR) to satisfy 
\cos{TB,TR)\ < 0.8. We retain 17.1% and 20.1% of lon- 
gitudinal and transverse signal, respectively, by requiring 
the aforementioned selection criteria. 

We calculate At = Az/P^c from the measured separa- 
tion Az between the Sj-ec and Btag vertices [40]. We de- 
termine the i?rec vertex from the two chargcd-pion tracks 
in its decay. The i?tag decay vertex is obtained by fitting 
the other tracks in the event, with constraints from the 
Brec momentum and the beam-spot location. The RMS 
resolution on At is 1.1 ps. We only use events that sat- 
isfy \At\ < 15 ps and for which the error on At {a At) is 
less than 2.5 ps. 

Signal events are identified using two kinematic vari- 
ables, the difference AE between the CM energy of the 
Brec, E^, and y/s/2, 

AE = E*s- v^/2, (14) 

and the beam-energy-substituted mass, 

ruES = y^(s/2 + P»-PB)V£^f-pL (15) 

where ^/s is the total CM energy. The Brec momentum 
Pb and four- momentum of the initial state {Ei,pi) are 
defined in the laboratory frame. We accept candidates 
that satisfy 5.25 < mss < 5.29 GeV/c^ and -0.12 < 
AE < 0.15 GeV. An asymmetric AE selection is used 
in order to reduce backgrounds from higher-multiplicity 
B decays. The resolution on niES is dominated by the 
knowledge of the energy of the e+ and e~ beams, and 
the resolution on AE is dominated by the reconstruction 
performance of the EMC. 

After the selection criteria mentioned above have been 
applied, the average number of candidates per event is 
approximately 1.5. In events with more than one can- 
didate, we select the candidate that minimizes the sum 
of {m-y-y — m„o)^ where m^o is the true 7r° mass [28]. In 
0.3% of events, the same tt" mesons are used by multiple 
B candidates. In such an event we randomly select the 
candidate to keep. 



IV. CONTINUUM BACKGROUND 
SUPPRESSION 

In addition to the cuts on cos 9i, R2, and 
I cos(TB,Ti?)| that directly remove background events, 

we use an artificial neural network in order to discrimi- 
nate between continuum background and signal events. 



For this purpose we combine the following eight variables 

into a single output, Af. 

• The coeflBcients, Lq,L2, split into sums over the 
ROE for neutral and charged particles; Lo^n,L'2,n 
and io,c, -^2,c- The coefficients are defined as Lk = 
'^Pj\ cos{ipj)\'^ , where fc = 0, 2, pj is the particle 
momentum and ipj is the angle of the particle direc- 
tion relative to the thrust axis of the B candidate. 
Both pj and ijjj arc defined in the CM frame. 

• I cos{B, Z)\, the absolute value of the cosine of the 
angle between the direction of the B and z axis in 
the CM frame. This variable follows a sine squared 
distribution for BB events, whereas it is almost 
uniform for qq. 

• I cos(TS, rj?)|. This variable, previously defined 

in Section IIIC, is strongly peaked at unity for 
qq events. BB events are more isotropic as the B 
mesons are produced close to the kinematic thresh- 
old. 

• I cos(T_B, Z) I , the absolute value of the cosine of the 
angle between the B thrust and the z axis. 

• The scalar sum of the transverse momenta pt in 
the ROE. This sum includes neutral and charged 
particles. 

The distributions of these input variables are shown in 
Fig. 4. Figure 5 shows Af for signal MC simulated events 
and continuum background samples (off-peak data) and 
the efficiency for signal and qq background as a function 
of cut on the minimum value of J\f. We require Af to be 
greater than —0.4. Note that later, we use this variable 
in the maximum likelihood fit described in Section VII. 

The samples used to train the neural net wore cor- 
rectly reconstructed MC simulated events and off-peak 
data. To avoid over-training, we used an independent 
sample of these data {i.e. distinct from the sample used 
for the training) to check the performance of the network. 
The training is stopped when the separation between the 
signal and continuum background is stable. 

V. MIS-RECONSTRUCTED SIGNAL AND 
SELECTION EFFICIENCY 

Mis-reconstructed signal candidates, referred to as self- 
cross-feed (SCF) signal, may pass the selection require- 
ments even if one or more of the pious assigned to the 
p+p" state belongs to the other B in the event. These 
SCF candidates constitute 50.7% (27.9%) of the accepted 
longitudinally (transversely) polarized signal. The ma- 
jority of SCF events have both charged pions from the 
p+p" final state, and unbiased CP information. These 
correct (right) track SCF events are denoted by RT SCF. 
There is a SCF component (13.8% of the signal) where 
at least one track in Srec is from the ROE. These wrong 
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I cos(TB, Z) I ^ P, (GeV/c) 



FIG. 4: The input variables used in training the neural network. The solid line represents the signal training sample (MC 
simulated events) and the dashed line represents the continuum background (ofT-peak data). The distributions shown are (in 
order left to right; top to bottom) I/o,n, Lo,c, L2,n, L2,c, cos{B, Z), cos(_B, TR), cos{TB, Z), and the sum of the pt in the ROE. 




Neural Network Output 



FIG. 5: The left plot shows the distribution of Af for off-peak data (dashed) and longitudinally polarized signal MC simulated 
events (solid). The right plot shows the signal efficiency (solid) and continuum background efficiency (dashed) distribution as 
a function of cut on the minimum value of TV. 



track (WT) events are used to determine the signal yield 
^signal and /i but have biased CP information, and are 
treated separately for the CP result. A systematic er- 
ror is assigned to the CP results from this type of signal 
event. The total selection efficiency for longitudinally 
(transversely) polarized signal is 7.7% (10.5%). 



VI. B MESON FLAVOR TAGGING 

A key ingredient in the measurement of time- 
dependent CP asymmetries is to determine whether at 
the time of decay of the Stag, corresponding to At = 0, 
the i?rec was a i?'^ or a 5°. This 'flavor tagging' is 
achieved with the analysis of the decay products of the 
recoiling B meson -Btag- The overwhelming majority of 
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B mesons decay to a final state that is flavor-specific, 
i.e. only accessible from either a B or a but not 
from both. The purpose of the flavor tagging algorithm 
is to determine the flavor of Stag with the highest pos- 
sible efficiency etag and lowest possible probability w of 
assigning a wrong flavor to Stag- It is not necessary to 
fully reconstruct i?tag in order to determine its flavor. 
In the presence of a finite mistag probability w, the CP 
asymmetry is reduced by a dilution factor l — 2w, so Eq. 
(7) becomes 

-lAtl/rgo f 

UiAt) = — 1±(1-2m;)x 

4tbo 



r]f[Ssm{AmdAt) - Ccos(AmdAt)]|. 

(16) 

The figure of merit for the performance of the tagging 
algorithm is the effective tagging efBciency 



Q = etag(l - 2wf 



(17) 



which is related to the statistical uncertainty a in the 
coefficients S and C through 



(J oc 



1 

7^' 



(18) 



We use a neural network based technique [11, 40] that iso- 
lates primary leptons, kaons and pions from B decays to 
final states containing D* mesons, and high momentum 
charged particles from B decays, to determine the flavor 
of the Stag- The output of this algorithm is divided into 
seven mutually-exclusive categories. These are (in or- 
der of decreasing signal purity) Lepton, Kaon I, Kaon II, 
Kaon-Pion, Pion, Other and Untagged. The performance 
of this algorithm is determined using fully reconstructed 
neutral B decays to D(*'~(7r+, p+, o^) and is summa- 
rized in Table II. The categories assigned correspond to 
events with leptons, kaons and pions in the decay prod- 
ucts of Stag- The Untagged category of events contain 
no flavor information and therefore carry no weight in 
the time-dependent analysis. In addition, these events 
are not considered useful for the branching fraction mea- 
surement since they are dominated by continuum back- 
ground. The total Q of this algorithm is 32.7 ± 0.7%. 



VII. LIKELIHOOD FIT METHOD 

On applying the selection criteria described above, we 
obtain a sample of 37424 events that enter the fit. These 
events are dominated by backgrounds from qq (81.4%) 
and BB (16.6%) events. The remaining 2% of events are 
considered as signal. We distinguish between the follow- 
ing components in the fit: 



TABLE II: Tagging efficiency e, average mistag fraction w, 
mistag fraction difference between and B° tagged events 



Aw, and Q for 



■ p events. 



Category 


ftag (%) 


to (%) 


Aij (%) 


Q (%) 


Lepton 


8.2 ± 0.1 


3.2 ± 0.5 


-0.2 ± 0.8 


7.2 ± 0.2 


Kaon 1 


11.3 ± 0.1 


3.7 ± 0.7 


1.1 ± 1.2 


9.7 ± 0.3 


Kaon II 


17.3 ± 0.2 


14.2 ± 0.7 


-0.9 ± 1.1 


8.8 ± 0.3 


Kaon-Pion 


13.4 ± 0.1 


20.8 ± 0.8 


0.5 ± 1.3 


4.6 ± 0.3 


Pion 


13.8 ± 0.2 


30.6 ± 0.8 


4.1 ± 1.3 


2.1 ± 0.2 


Other 


9.4 ± 0.1 


40.1 ± 1.0 


2.3 ± 1.5 


0.4 ± 0.1 


Untagged 


26.8 ± 0.2 50.0 ± 0.0 




0.0 ± 0.0 


Total 


32.7 ± 0.7 



(ii) SCF signal, split into RT and WT parts, 

(iii) charm B"^ backgrounds (6 — > c), 

(iv) charm B° backgrounds (6 — > c), 

(v) charmless B^ backgrounds, 

(vi) charmless backgrounds, 

(vii) continuum background. 

The dominant B backgrounds come from components 

(iii) and (iv). The signal, continuum and B background 
models are described in Sections Vll A, Vll B, and Vll C, 
respectively. 

We use an unbinned, extended maximum likelihood 
(ML) fit to extract A^signai, /l, -Siong, and Ciong- The 
likelihood function incorporates the following previously 
defined eight discriminating variables to distinguish sig- 
nal from the backgrounds: TOes, AE, At, and the 
and cos Ok values of the two p mesons, where; fc = 1, 2. For 
each of the aforementioned components j we construct a 
probability density function (PDF) that is the product 
of one-dimensional PDFs for each of the variables. The 
PDFs do not account for all possible correlations among 
the discriminating variables and we account for possible 
biases as a result of neglecting these correlations as dis- 
cussed in Section Vlll. For each event i, the PDFs can 
be written as 

= V^{rr^Es)'Pj{^E')-Pj{^n'PA^') x 

Vj (mi )Vi {m\)Vi (cos e\ )Vi (cos 4) . (19) 

We determine the PDFs for signal and each of the B 
background components from MC distributions for each 
discriminating variable. The PDFs for the continuum 
background are determined from on-peak and off-peak 
data. The likelihood function is 



C = 



exp(- Ej %) 

m 



N 



(20) 



(i) correctly reconstructed signal, 



where rij are the event yields of hypotheses j (determined 
from the ML fit) and N is the observed number of events 
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in the sample. The normahzed exponent takes into ac- 
count Poisson fluctuations in tlie total number of events. 
We simultaneously fit events in the mutually-exclusive 
flavor tag categories listed in Table 11, excluding events 
in the Untagged category. There are six continuum back- 
ground yields, one for each flavor tag category, and we 
use a single yield for the signal and each of the B back- 
grounds, while accounting for the relative fractions of 
events expected in each of the flavor tag categories. The 
fit results are obtained by minimizing the value of — ln£. 



A. Signal model 



The signal has contributions from true and SCF longi- 
tudinally (long) and transversely (tran) polarized events. 
In addition to this, the longitudinally polarized SCF sig- 
nal PDF is further sub-divided into the two categories of 
mis-reconstnictcd signal: RT and WT SCF signal, and 
all signal PDFs are sub-divided according to the i?tag 
flavor tag category. The total signal PDF is given by 



signal — JL\[i^ ^RTSCF ^ WT SCfJ '^true 



I f long -Tjlong . j^long T->long 
"•"/RTSCF^'^RTSCF /Vi^TSCF^WTSCF 



+(1 - urn - frsF]Ku: + /rcWc?)- 



(21) 



where /^^^^^ is the fraction of longitudinally po- 
larized RT(WT) SCF signal and f^}^§ is the frac- 
tion of transversely polarized SCF signal. The PDFs 
V^^^-^j,^ gpp and T^gcF ^'^'^ defined accordingly. In order 
to extract the observable /i, from the fit, we account for 
the diflferent reconstruction efficiencies for longitudinally 
and transversely polarized signal using 



fL = 



£ohs 
JL 



(22) 



where is the fraction of signal from longitudinally 
polarized events in the data, and R is the ratio of longi- 
tudinal to transverse polarized signal efficiencies. 

All of the SCF PDFs are further sub-divided into 
parts that contain correctly reconstructed and mis- 
reconstructed p mesons. Four sub-divisions are used, cor- 
responding to both mesons being correctly reconstructed 
{TT), one being correctly reconstructed and the other 
being mis-reconstructed {TF or FT), and both mesons 
being mis-reconstructed (FF). Here the first (second) p 
meson is positively (negatively) charged. The SCF PDFs 



have the form 



V = V{M)V{M)[ 



frrVrT [mES , ^E, cos 6*1 , toi , cos 92,m2) + 
fTpVTF {mES , AS, cos 6*1 , TOi , cos 62,m2) + 
IftVet (mES, ^E, cos 9i,mi, cos 6*2, m2) -I- 
fpFVpFimES, ^E, cos ^1, mi, cos ^2, ^2)], 

(23) 

where T and F are correctly and mis-reconstructed p 
mesons, respectively, /tt, /tf, ./ft, Iff are the frac- 
tions of each type of mis-reconstructed event, and the 
PDFs for niES, ^E, and the helicity and mass for each 
reconstructed p meson are the products of one dimen- 
sional PDFs. The signal PDFs for TV and At have dif- 
ferent shapes according to the Stag flavor tag category 
assigned to an event. 

The longitudinally polarized true and RT SCF signal 
niES distribution are described by a Gaussian with a 
power law tail [41-43] which takes the following form 

\ ( {m'-mof\ 
f{m) oc exp ^ ) , m>mo-^a, 



oc 



2a2 

n" exp(-^V2) 



1) 



n — 



- , m <mo — ^a, (24) 



where toq is the mean, a is the width of the Gaussian 
part, and the parameters ^ and n describe the exponen- 
tial tail. The remaining PDFs for the signal niES distri- 
butions are the sum of a Gaussian and a Gaussian with 
an exponential tail. The signal AE distributions for the 
longitudinally polarized SCF events arc described by first 
and second order polynomials, with the exception of TT 
events, which are modeled with the sum of a second or- 
der polynomial and a Gaussian. The other IS.E PDFs are 
described by the sum of Gaussian and a Gaussian with 
an exponential tail. The M distributions are modeled us- 
ing smoothed histograms of MC simulated events. The 
signal Ai distribution is described by Eq. (16) convolved 
with a triple Gaussian resolution function given by 

7?.sig(Af, (7Ai) = /coroG (At, ^coretTAt, tTcore^At) + 
/tailG (At, /itailCAt, (^tai\(^ At) + 
/outlierG (At, /^outlier, foutlier) , (25) 

where G is a Gaussian with means pbi and width <Ji for 
I = core, tail, and outlier. The parameters Utaii, o'outiier 
and ^outlier are set to 3.0 ps, 8.0 ps and 0.0 ps, re- 
spectively. The remaining parameters of the resolution 
function are determined from signal MC simulated events 
scaled by the differences between data and MC simulated 
events for large samples of exclusively reconstructed B 
decays to £'^*^~(7r"'", p"*", a^) final states described in Sec- 
tion VI. The values of the mean and width of the core 
Gaussian contribution to the resolution function depend 
on the flavor tagging category of an event. There are 
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three signal At PDFs, one for the true and RT SCF lon- 
gitudinally polarized signal, one for the WT SCF longitu- 
dinally polarized signal and one for the transversely po- 
larized signal. Each of these PDFs has distinct values of 
S and C as described in Sections I and V. The WT SCF 
longitudinally polarized signal has a different resolution 
function with respec;t to the rest of the signal. The cos 6i 
distribution for true signal events is described by the ex- 
pected cosine squared or sine squared distribution mul- 
tiplied by a polynomial acceptance function. The cos^j 
PDFs for T p mesons in longitudinally (transversely) po- 
larized SCF signal are described by sixth order polyno- 
mials (smoothed histograms of MC simulated events). 
The cos 9i PDFs for mis-reconstructed p mesons are de- 
scribed by smoothed histograms of MC simulated events. 
The longitudinally polarized signal p mass PDFs are de- 
scribed using relativistic Breit-Wigners and third order 
polynomials for T and F p contributions, respectively. 
The transversely polarized signal SCF p mass PDF dis- 
tributions are described using smoothed histograms of 
MC simulated events. 

In addition to A^signab /l) 'S'long) and Ciong, we deter- 
mine the mean of the core Gaussian part of the m es and 
AE distributions from the fit. We obtain parameters 
consistent with expectations from MC simulated events. 



B. Continuum background model 

The continuum background PDFs are also sub-divided 
according to the mis-reconstruction of p mesons in the 
final state and have the form 



r = 



V{mEs)V{AE)V{Af)V{At) 
/tt'Ptt (cos 



9i,mi,cos 62,1712) 
6i,mi,cos6: 
61, mi, cos 9: 

fp pPp p (cos Oi,mi, cos 9[ 



fj.pPj,p(^cos 9i, mi, cos 92,1712) + 
IftVpt (cos 9i,mi, cos ^2 , ^2) + 

^2, "12) 



(26) 



where helicity and mass PDFs for each reconstructed p 
meson are the products of one dimensional PDFs. 

The continuum distribution for mps is described by a 
phasc-space-motivated distribution [44] with the follow- 
ing form 

f{x) (X xVl-a;2 • exp[e (1 - x")] ■ 9{mEs), (27) 

where 9{mEs) = 1 for mps < v^/S and 9{mEs) = 
for mES > \fs/'^ and x = 2mEsl \fs- The AE and M 
shapes are modeled with third and fourth order polyno- 
mials, respectively. The parameters of the mEs, AE and 
Af shapes are allowed to vary in the fit to the on-peak 
data. The continuum p mass distribution is described us- 
ing a relativistic Breit-Wigner and a third order polyno- 
mial PDF for T and F p contributions, respectively. The 



cos 9i distribution is described by a third order polyno- 
mial. The continuum At distribution has a prompt life- 
time component convolved with a triple Gaussian resolu- 
tion function. The parameters of the p mass and helicity 
distributions are obtained from a fit to the off-peak data, 
and the remaining parameters are determined in the fit. 



C. B background model 

1. Charm B backgrounds 

Combinatorial events from b ^ c transitions are 
the dominant B backgrounds. These components have 
shapes similar to continuum and do not peak in the signal 
region for the discriminating variables. The functional 
form used for the PDFs of these background components 
is given by Eq. (26). We parameterize the m.Es and 
AE distributions of these backgrounds using the phase- 
space-motivated distribution of Eq. (27), and a third 
order polynomial, respectively. The p mass distribu- 
tion is described using a relativistic Breit-Wigner and a 
third order polynomial PDF for T and F p contributions, 
respectively. The remaining PDFs are described using 
smoothed histograms of MC simulated events. Each of 
the PDF parameters for these backgroimds are deter- 
mined from samples of MC simulated events, and the 
yields of these components are determined in the fit. 
When studying systematic uncertainties we modify the 
PDF used for At so that it has a form similar to the sig- 
nal which uses the resolution function of Eq. (25). An 
effective lifetime that is smaller than Tgo is used to ac- 
count for mis-reconstruction of these events and the finite 
charm meson lifetime. The value of this parameter is ob- 
tained by fitting MC simulated events for this category 
of events. 



2. Charmless B backgrounds 

Some of the charmless B backgrounds have PDFs sim- 
ilar to the signal for one or more of the discriminating 
variables, so it is important to correctly account for such 
events in the fit. We consider the 20 exclusive and 2 in- 
clusive components of this type of B background listed 
in Table III. If a charmless B background decay con- 
tributes an event yield more than 1% of the expected 
signal yield, we model that mode exclusively. The re- 
maining 140 charmless background components consid- 
ered were combined with the appropriate weightings to 
form the neutral and charged inclusive charmless B back- 
ground components. Where possible we use branching 
fractions from existing measurements. Where measure- 
ments are not available, we have either tried to use SU (2) 
and SU (3) flavor symmetries to relate an unmeasured 
decay to a measured one, or where this is not possible, 
we have assigned a branching fraction of 10 xlO~^ to a 
decay mode. An uncertainty of 100% is assumed on all 
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TABLE III: The components of charmless B backgrounds 
considered, along with the branching fraction (B) and number 
of events (Ntg) expected in the selected data sample, where 
f indicates a longitudinally polarized final state, and NR de- 
notes a non-resonant final state. Where appropriate, branch- 
ing fractions of these decays are quoted, including the branch- 
ing fraction of sub-decay modes and measured or expected 
values of Jl- 



Decay Mode 


BxlO"^ 


Nbg 


B+ 




K+-K+n-[45, 46] 


00 ± O 


11 rt i 


B+ 






10 ± 10 


15 ± 15 


B+ 




7r'V+ [47-49] 


11 1 o 

ii ± Z 


/o ± iU 


B+ 




p+p°[21, 22] t 


17 lb 3 


71 ± 12 


B+ 








O -4- Q 


B+ 




K|(1770)p 


10 di 10 


9 ± 9 


B+ 




A',>[45, 46, 49] 


1 1 _i_ 1 
11 ± 1 


n _i_ 1 
9 ± 1 


R+ 




0,^ TV 


on -1- on 


firi -1- (if\ 
DU ± OU 


B+ 






10 ± 10 


12 ± 12 


B+ 




a?(p+^-)p+t 


8 ± 8 


13 ± 13 


B+ 




a?7r+ 


20 ± 20 


49 ± 49 


B+ 




charmless 




59 ± 15 


B° 




7r+7r-7r° (NR) 


30 ± 3 


23 ± 23 


B° 




p±7r=F[49, 50] 


24 ± 3 


42 ± 4 


B° 






23 ± 23 


15 ± 15 


B° 




af (p%±)7r=F[51] 


20 ± 2 


7 ± 1 


B° 




p=F7r±7r« 


10 ± 10 


45 ± 45 


B° 




Jf2*(1770)p 


10 ± 10 


8 ± 8 


B° 




a±(p±^»)pT [52] t 


16 ± 16 


43 ± 43 


B° 




a±(pV)pT [52] t 


16 ± 15 


9 ± 8 


B" 




a±(p±7r°)7rT[51] 


40 ± 4 


102 ± 9 


B° 




charmless 




88 ± 22 



extrapolated branching fractions. Wc assume that the ai 
meson decays into a three pion final state via pn. The 
decay B'^ — > a?7r° is penguin dominated and is expected 
to have a small branching fraction, as is the case for the 
penguin dominated decays B° p°7r° and P^P°- 
As a result, we have neglected possible contributions from 
this potential source of background. Although we don't 
see evidence for B° p^tt^tt^, we do consider this 
channel as a potential source of background. Contribu- 
tion from Att is expected to be even smaller and 
is neglected. When considering the systematic uncer- 
tainty arising from interference between signal and other 
7r"'"7r~7r°7r° final states, we assume that the non-resonant 
7r+7r~7r'^7r° final state has the same branching fraction 
as that of B" p'^tt^tt", where B(B" p'^Tr^Tr") is 
calculated from the yield determined in the B^ P'^P~ 
nominal fit. This is discussed further in the Appendix. 

The functional form used for the PDFs of these charm- 
less i? background components is given by Eq. (26). The 
AE distributions are described by third order polynomi- 
als, except for non-resonant B^ p^tt^tt'^ which uses 
smoothed histograms of MC simulated events. The mp 
distributions for true p mesons are parameterized using a 
rclativistic Breit-Wigner, and the fake nip (combinatorial 
7r^7r°) distribution is described using third order polyno- 



mials. The remaining background shapes are described 
using smoothed histograms of MC simulated events. The 
yield for i?" p^n^n^ decays is allowed to vary in the 
fit. All other charmless background yields are fixed to ex- 
pectations. This constraint is relaxed when studying pos- 
sible sources of systematic uncertainties. When studying 
systematic uncertainties from possible CP violation in 
the B background, we modify the PDF used for At so 
that it has a form similar to the signal which uses the 
signal resolution function in Eq. (25). 



VIII. LIKELIHOOD FIT VALIDATION 

Before applying the fitting procedure to the data, we 
subject it to various tests. The aim of these tests are to 
verify that one can extract the signal observables A^signai, 
/l, iSiong, and Ciong in a controlled way. Consistency 
of the likelihood fit is checked with ensembles of exper- 
iments simulated from the PDFs. The event yields are 
generated according to a Poisson distribution with mean 
rij for each category in these ensembles. In each of these 
tests we verify that the values generated for the signal 
observables are reproduced with the expected resolution. 
The distribution of — In £ for the ensemble of experi- 
ments in comparison to that obtained when fitting the 
data provides an additional, but limited, cross-check of 
the consistency. 

The PDFs used in the likelihood do not account for all 
possible correlations among the discriminating variables. 
We account for possible biases as a result of neglecting 
these correlations by fitting ensembles of experiments ob- 
tained from samples of signal and the charmless B back- 
ground MC simulated events combined with charm back- 
grounds and qq background events simulated from the 
PDFs. The MC simulated events used in these ensem- 
bles have these correlations modeled. We find a positive 
bias of 38.5 ± 5.3 events on iVgignai, and a negative bias of 
0.016 ± 0.002 on /l, and we do not observe a significant 
bias on 5iong and Ciong- 

As continuum events are the dominant background, 
we apply the fitting procedure to the off-peak data (after 
correcting for the difference in ^/s) to confirm that we do 
not find a fake signal in this control sample of events. We 
fit 790 off-peak events and extract signal and continuum 
yields of 8 ± 7 and 782 ± 28 events, respectively. 

A blind analysis technique has been used for the ex- 
traction of ^long and Ciong, where the actual values of 
these observables have been hidden by an offset. With 
the values of S'long and Ciong hidden we perform the fol- 
lowing fit cross-checks. We first verify that the uiic;er- 
tainties on the signal observables, and the value of — In £ 
obtained from the fit to data, are compatible with the en- 
sembles of simulated experiments described above. We 
then validate the stability of our results by introducing 
a variety of modifications to the fitting procedure. In 
addition to the nominal set of variables determined from 
data, we allow the physics parameters tbo and Anid to 
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vary in turn. The signal obscrvablcs do not change signif- 
icantly when doing this, and the results obtained for t^o 
and Arud are 1.72±0.16 ps and 0.36±0.22 7i/ ps, respec- 
tively, consistent with the reported world average [28] . 

The Lepton and Kaon I tagged events have the highest 
signal purity, and lowest w probability, and dominate 
our CP results. We check that the results obtained from 
fitting only these categories of events are compatible with 
our nominal blind results. When doing this, we observe 
shifts of —0.08 (—0.05) on S\ong (Ciong) relative to the 
nominal result. 

Given that there are a number of B backgrounds that 
contribute to the data that have yields fixed to expecta- 
tions, we validate this assumption by allowing each fixed 
B background yield to vary in the fit. We obtain back- 
ground yields consistent with our expectations, and ob- 
serve the shifts on signal parameters listed in Table IV. A 
systematic uncertainty corresponding to the largest ob- 
served deviation is assigned on our results. 

Once these checks have been completed, the fitting pro- 
cedure is frozen and we extract the values of S'long and 
C'long by removing the offset. 



TABLE IV: The shifts on signal yield ((5iV(signal)), fraction 
of longitudinally polarized events (<5/l), Siong (SSiong), and 
Ciong (5Ciong) obtained when floating the yield of each fixed B 
background in turn. The symbol f indicates a longitudinally 
polarized final state, and NR denotes a non-resonant final 
state. The yields of decay modes not listed in this table are 
allowed to vary in the nominal fit. 



Decay Mode 


(5A/'(signal) 




*5 S'long 


(5 Ciong 


B+ K+n+n- 


3 


-0.002 


-0.002 


-0.003 


B+ ^ w+7r\° 


11 


0.001 


0.002 


-0.001 


B+ ^ 7T°p+ 


9 


0.000 


0.001 


-0.002 




-8 


-0.001 


-0.002 


0.000 


B+ p+TT+n- 


1 


0.000 


-0.001 


-0.001 


B+ K2{1770)p 


-5 


0.001 


-0.001 


-0.000 


B+ -> KStt 


-13 


-0.001 


-0.002 


0.001 


B+ ^ a+7r° 


17 


0.003 


0.006 


-0.001 


B+ ^ p+7v\° 


-2 


0.000 


0.001 


-0.002 


B+ ^ a;(p+7r-)p+t 


-34 


-0.006 


-0.028 


-0.003 


B+ ^ a?7r+ 


5 


0.001 


0.004 


-0.001 


B"*" — » charmless 


-3 


-0.001 


-0.002 


0.000 


B° ^ 7r+7r-7r° (NR) 


-4 


0.000 


0.000 


0.000 


B° ^ p±7r=F 


-5 


0.000 


-0.001 


0.000 




-25 


-0.003 


-0.009 


-0.001 


B°^a±(p°7r±)7r=F 


1 


-0.002 


-0.001 


-0.002 


B° K^{1770)p 


-31 


-0.003 


-0.009 


-0.003 


B° ^ a±(p±7r°)pTt 


-25 


-0.002 


-0.025 


-0.004 


B°-ahpV*)p^t 


11 


0.000 


0.006 


0.000 


B°^a±(p±7r")7rT 


10 


0.002 


0.003 


-0.002 


B° — » charmless 


1 


0.000 


0.000 


-0.001 



IX. RESULTS 

A. Fit results 

From the ML fit described above, we obtain the fol- 
lowing results 

A''(signal) = 729 ± 60(stat), 

/i = 0.992 ± 0.024(stat), 

5iong = -0.17±0.20(stat), 

Ciong = 0.01±0.15(stat), 

after correction for a -1-38.5 ± 5.3 event fit bias (see Sec- 
tion VIII), a —76 event bias from SCF on the signal yield 
(see Section IX B) and a correction for a —0.016 ± 0.002 
fit bias on f^. The correlation between S'long and C'long 
is —0.035. We discuss the origin of these fit biases in 
Section IX B. We calculate the branching fraction of this 
decay using B = N signal / {^N^^-j^)^ where e is the efficiency 
for signal corresponding to the observed /l, and Nb pairs 
is the number of BB pairs analyzed. We obtain 

6(S"^p+p^) = (25.5±2.1(stat))xl0-^ 

The i?° p'^TT^TT^ background yield obtained from the 
fit is 9.2 ± 53.6 events. Figure 6 shows distributions of 
niES, AE, cos 9i and m^±^o for the Lepton and Kaon 
I tagged events with a loose requirement on J\f. Rela- 
tive to the total number of events in the fit, the plot of 
iTiEs contains 15.6% of the signal and 1.1% of the to- 
tal background. For the other plots there is an added 
constraint that tues > 5.27GeV/c^; these requirements 
retain 13.9% of the signal and 0.4% of the total back- 
ground. Figure 7 shows the At distribution for 5° and 
B^ tagged events, as well as the time-dependent decay- 
rate asymmetry of Eq. (10). Here we apply the same 
selection criteria as in Fig. 6(b)-(d). 

B. Systematic uncertainty studies 

Table V lists the possible sources of systematic uncer- 
tainties on the values of the A^signah /l, S'long, and Ciong 
that have been studied, and are described in the follow- 
ing. 

• The uncertainty from PDF parameterisation is ob- 
tained by varying PDF shape parameters by ±lcr, 
in turn. The deviations obtained are added in 
quadrature to give the quoted uncertainty. 

• The systematic uncertainty from the fraction of 
SCF events is obtained from the difference between 
the nominal fit result and that obtained when fit- 
ting the data and removing the SCF from the fit 
model. In the case of the signal yield we correct 
for half of the difference observed, and use 100% of 
the correction as an uncertainty. The uncertainty 
on the other signal observables comes from the dif- 
ference observed between the two sets of fit results. 
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TABLE V; Summary of additive systematic uncertainty contributions on the signal yield A'signai, /l, 5iong and Ci, 



Contribution 


0"(Afsignal) 




O"(5'long) 


""(Clong) 


PDF parameterisation 


+ 23 

-41 


+0.019 
-0.005 


+0.02 
-0.04 


0.03 


SCF fraction 


76 


0.003 


0.00 


0.02 


Control samples calibration 


16 


0.006 


0.01 


0.01 


rriES endpoint 


12 


0.001 


0.00 


0.01 


B background normalization 


+ 16 
-20 


+0.005 
-0.002 


0.01 


0.00 


floating B backgrounds 


34 


0.006 


0.03 


0.00 


B background CP asymmetry 


2 


+0.001 
-0.000 


0.00 


+0.02 
-0.01 


Tb 


+0 

-1 


+0.001 
-0.000 


0.00 


0.00 


ArUd 


+0 

-1 


+0.001 
-0.000 


0.00 


0.00 


tagging and dilution 


+2 
-16 


+0.010 
-0.001 


0.00 


0.01 


transverse polarization CP asymmetry 


+0 
-9 


+0.006 
-0.000 


0.01 


0.01 


Wrong track SCF CP asymmetry 


+0 
-3 


+0.001 
-0.000 


0.01 


0.01 


DCSD decays 






0.01 


0.04 


Interference 


18 


0.000 


0.01 


0.01 


Fit Bias 


19 


0.008 


0.02 


0.02 


SVT Alignment 






0.01 


0.01 


Total 


+94 
-102 


+0.03 
-0.01 


+0.05 
-0.06 


0.06 




5.25 5.26 5.27 ,5.28 5.29 
rrigg (GeV/c") 




-0.5 0.5 
cos(6,) 




30 
20 
10 

30 
20 
10 
1 

0.5 


-0.5 
-1, 











(a) 










[ (b) 




r 
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m,,,. (GeV/c-) 



-2 
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FIG. 6: The distributions for the highest purity tagged events 
for the variables (a) rriES, (b) Ai?, (c) cosine of the p helic- 
ity angle and (d) m^±^o. The dashed lines are the sum of 
backgrounds, and the solid lines are the full PDF. 



FIG. 7: The At distribution for a sample of events enriched in 
signal for (a) B'' and (b) B^ tagged events. The dashed lines 
are the sum of backgrounds, and the solid lines are the sum of 
signal and backgrounds. The time-dependent CP asymmetry 
(see text) is shown in (c), where the curve is the measured 
asymmetry. 



• The kinematic endpoint position in is ex- 

tracted from the fit. Changes in beam energy in 
the data can affect the endpoint position. To ac- 



count for possible effects of this, we vary the kine- 
matic endpoint position in mES by ±0.45 MeV/c^ 
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which is determined from control samples of fully 

reconstructed i?° decays. 

The uncertainty from the niES and AE widths is 
obtained from the observed shifts relative to our 
nominal result, when allowing these parameters to 
vary independently in the fit to data. 

We vary the B background normalization within 
expectations for each background in turn. The de- 
viations obtained arc added in quadrature to give 
the quoted uncertainty from this source. 

As the branching fractions of some of the B back- 
grounds are not well known, we assign an additional 
uncertainty coming from the maximum shifts ob- 
tained when allowing each of the fixed backgrounds 
to vary in turn in the fit to data. 

Additional uncertainties on the CP results come 

from possible CP violation in the B backgroimd. 
We replace the At PDFs used by each of the B 
backgrounds in turn by one resembling the signal. 
Charged B backgrounds can have non-zero values 
of Acp, and neutral B backgrounds can have non- 
zero values of S and C. We use existing experimen- 
tal constraints where possible, otherwise we allow 
for a CP asymmetry up to 10% in B decays to final 
states with charm, and up to 50% in B decays to 
charmless final states. 



• As the PDFs used in the ML fit do not account for 
all of the correlations between discriminating vari- 
ables used in the fit, the results have a small bias. 
We calculate the fit bias on the signal observables 
as described in Section Vlll. The observed bias on 
-^signal, and /l is corrected, and 100% of the cor- 
rection is assigned as a systematic uncertainty. Wc 
do not observe a significant bias on Siong and Ciong, 
and conservatively assign a systematic uncertainty 
from fit bias based on the statistical precision of 
this test. 

• Small imperfections in the knowledge of the geom- 
etry of the SVT over time can affect the measiirc- 
ment of Siong and Ciong- We vary the alignment 
according to the results obtained from the study of 
e~^e~ e^e^, events in order to estimate 
the magnitude of this systematic error on our CP 
results. 



The branching fraction has multiplicative systematic 
uncertainties from the reconstruction of 7r° mesons in 
the detector (6%), uncertainties in the reconstruction of 
charged particles (0.8%), and the discrimination of tt^ 
from other types of charged particles (1%). In addition 
to these uncertainties, there is a 1.1% uncertainty on the 
number of BB pairs in the data sample. The statistical 
uncertainties arising from the MC samples used in this 
analysis are negligible. 



The physics parameters tbo = 1.530 ±0.009 ps and 
Arud = 0.507 ± 0.005 h/ ps [28] are varied within 
the quoted uncertainty. 

The tagging and mistag fractions for signal and 
the B backgrounds are corrected for data/MC dif- 
ferences observed in samples of fully reconstructed 
hadronic B decays. Each of the tagging and mistag 
parameters is varied in turn by the uncertainty 
from the correction. The deviations obtained are 
added in quadrature to give the quoted uncertainty 
from this source. 

Allowing for possible CP violation in the transverse 
polarization, and in the WT longitudinally polar- 
ized signal SCF events results in additional uncer- 
tainties on signal observables. We vary S and C by 
±0.5 (±1.0) for the transverse polarization (WT 
SCF). 

Possible CP violation from interference in doubly 
Cabibbo-suppressed decays (DCSD) on the tag side 
of the event [53] contribute to systematic uncertain- 
ties on iSiong and Ciong- 

We estimate the systematic error on our results 
coming from neglecting the interference between 
B^ — > p'^p~ and other 47r final states: B — > ai7r, 
/Otttt" and tt-k-k^tt^ . Details of this calculation are 
given in the Appendix. 



X. CONSTRAINTS ON THE UNITARITY 
TRIANGLE 

A. The i5l7(2) isospin analysis 

In SU (2) isospin symmetry, the amplitudes of B decays 
to pp final states [20] are related by: 



72 
\/2 



00 



^+0 - A 



-r-O -rOO 

■ A - A , 



(28) 
(29) 



for the longitudinal polarization and each of the CP 

eigenstates of the transverse polarization, where A^^ 

(A*"*) are the amplitudes of B{B) decays to the final state 
with charge ij, where i = +,—,0 and j = — ,0. These 
two relations correspond to triangles in a complex plane 
as shown in Fig. 8. In the usual phase convention [20], 
the amplitudes A are rotated to A in order to align the 
base of the triangles. After aligning A+° and A , the 
phase difference between A'^ and A'^ is 26a. There 
are two sources of ambiguities on the measurement of a. 
There is a four-fold ambiguity coming from the orienta- 
tion of the isospin triangles, and a two-fold ambiguity 
from the trigonometric conversion in Eq. (12). In order 
to measure a one must measure the branching fractions 
and charge asymmetries of iJ decays to p'^p~ , p^p^, P^P°- 
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FIG. 8: The isospin triangle for B —> pp decays. 

There are several assumptions implicitly used in the 
isospin-based direct measurements of a: 

• This approach only considers tree and gluonic pen- 
guin contributions. Possible contributions from 
electroweak penguins (EWP) are neglected as they 
do not obey SU{2) isospin symmetry. These EWPs 
have the same topology as the gluonic penguin di- 
agram in Fig. 2, with the gluon replaced by 7 or 
Z° bosons. In the absence of EWP contributions 
1^+0 1 = \A~°\, and Acp is zero for B+ ^ p+ p° . 
Several groups have estimated the correction due to 
the SU{2) breaking effect of EWP contributions to 
be 1.5 — 2.0 degrees [54, 55]. These estimates con- 
sider contributions from the two EWP operators 
assumed to be dominant in the effective Hamilto- 
nian. 

• The possible effect on the isospin analysis from p'^ — 
uj mixing [55] is neglected. 

• Other SU{2) symmetry breaking effects are ne- 
glected. Estimates of the magnitude of these effects 
are much less than the current experimental preci- 
sion [56, 57]. Possible isospin 5/2 amplitudes also 
break the SU{2) triangle construction [58]. 

• The isospin analysis outlined above neglects possi- 
ble 1=1 amplitudes [59]. The presence of 1=1 am- 
plitudes in B° P^P^ can be tested by measuring 
"Siong and Ciong for different ranges of the invariant 
TT^TT^ mass. 

We constrain the CKM angle a and penguin contribu- 
tion Sa from an isospin analysis of B pp decays. The 
inputs to the isospin analysis are the amplitudes of the 
CP-even longitudinal polarization of the pp final state, 
as well as the measured values of 5iong and Ciong for 
B'^ p~^ p~ ■ We use the following numerical inputs in 
the isospin analysis: 

• The measurements of B{B'^ — > p'^p^), f^, S'long, 
and Ciong presented here. 

• The measurements of B and /l for P'^p'^ 
from Ref. [21]. 

• The measurement of P^ P^) from Ref. [23]. 




a (degrees) 



FIG. 9: Confidence level on a obtained from the isospin anal- 
ysis with the statistical method described in [14] . The dashed 
lines correspond to the 68% (top) and 90% (bottom) CL in- 
tervals. 



To interpret our results in terms of a constraint on 
a from the isospin relations, we construct a that in- 
cludes the measured quantities expressed as the lengths 
of the sides of the isospin triangles and we determine the 
minimum Xo- We have adopted a simulated-experiment 
technique to compute the confidence level (CL) on a; our 
method is similar to the approach proposed in Ref. [60]. 
For each value of a, scanned between and 180°, we 
determine the difference Axdata('^) between the min- 
imum of x^{a) and Xo- We then generate MC experi- 
ments around the central values obtained from the fit to 
data with the given value of a and we apply the same 
procedure. The fraction of these experiments in which 
Ax^(^(a) is smaller than IS.Xt)Ata{'^) interpreted as 
the CL on a. Figure 9 shows 1 — CL for a obtained from 
this method. 

It is possible to obtain a prediction of a from indi- 
rect constraints by combining measurements of the CKM 
matrix elements |Ks|, |Kd|, |Kb|, and \Vch\, CP vi- 
olation in mixing from neutral kaons, B-B mixing in 
Bd and Bg mesons and the measurement of sin 2/3 from 
b — > CCS decays. The indirect constraint on a from the 
UTfit [15] and CKMfitter [14] groups are (92.9 ± 5.7)° 
and (100l7;3)° respectively. Recent calculations using 
QCD Factorization predict a = (85.6+^^:3)° [61]. Select- 
ing the solution closest to the CKM combined fit average 
we find a = [73.1, 117.0]° at 68% CL, where the error is 
dominated by \5a\ which is 18° at 68% CL. The con- 
straint obtained on a is worse than the previous result in 
Ref. [26] for the following reasons: (i) the central value 
of the branching fraction for B^ P^ P^ [21] is smaller 
than the value previously used [22, 62], and (ii) the re- 
sults of the latest search for p^ p^ show evidence for 
a signal [23] , while the previous constraint on a used the 
central value corresponding to the upper limit available 
at that time [63] . Both of these factors lead to an increase 
in the penguin contribution to the total uncertainty on 
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a when using an isospin analysis. 

A future extension to the SU{2) isospin method to 
measure a described here will be possible when there are 
sufficient statistics to perform a time-dependent analysis 
of p^p^. The parameters S'long and Ciong measured 

in —f p^p^ can be incorporated into the isospin anal- 
ysis, allowing one to over-constrain the isospin triangles 
of Eq. (28) and Eq. (29). 

B. Flavor SU{3) analysis 

There has been progress on constraining the pen- 
guin contribution in — > p^ p^ using a flavor SU{3) 



based approach and experimental constraints from — > 
K*°p+ [24]. The amphtudc for i?° p+ p~ can be writ- 
ten as 



AiB"^ ^ p+p-) ^ Te''' + Pe'^-^" , (30) 



where T and P are the magnitudes of tree and penguin 
amplitudes, 7 is the phase of Vub, and (5tp is the strong 
phase difference between T and P. The approach of 
Ref. [24] relates the penguin amplitude in B'^ P^P^ to 
the amplitude of the penguin decay B^ —f K*'^p^ giving 
three relations: 



^long 

i4d|/p \'rL(B±^i^*V) 



2r sin^TP sin(/3 



1 — 2r cos Stp cos(/3 + a) +r^' 

sin 2a + 2r cos <5tp sin(/3 — a) — sin 2(3 



1 — 2r cos Stp cos(/3 



VcsIJk'J Ti^iBO^p+p-) l-2rcos(STPCOs(/3 + a)+r2' 



(31) 

(32) 
(33) 



with three unknowns, a, r, and (Stp- The parameter 
r = \P/T\, and (3 is the phase of Vtd- We use the 
value of (3 obtained from b ccs decays [7, 8]. The 
CP averaged decay rates of the longitudinal components 
of B^ K*°p+ and B° p+ p- (El) are related by 
the squared ratio of CKM matrix elements and decay 
constants (/i, where p or K*) as shown in Eq. (33). The 
factor F accounts for additional sources of SU{i) break- 
ing not described by the decay constants. 
The assumptions used in this approach are: 

• The amplitude in the penguin dominated decay 
5+ — ^ K*^p^ is related to the penguin amplitude 
in B° p+p-, 

• The dominant SU{3) breaking correction ac- 
counted for by F is the neglect of annihilation dia- 
grams in the B^ — > K*^p'^ decay. We use the same 
value of F as Ref. [24]; F = 0.9± 0.6. 

The result is a constraint on a with theoretical uncer- 
tainty from penguin contributions of (-1-2, —5)° which 
represents a stronger constraint on the unitarity triangle 
than provided by the isospin analysis constraint. Fig- 
ure 10 shows 1 — CL for a obtained with this method. Se- 
lecting the solution closest to the CKM combined fit av- 
erage [14, 15] we find a = [83.3, 105.8]° at 68% CL, where 
the error is dominated by the experimental uncertainty 
on S'long and Ciong- The strong phase difference Stp is 
only weakly constrained, and r = 0.10jiQ Q4. The QCD 



factorization approach in Refs. [64-66] predicts that Stp 
is small, and so we require \6tp\ < 90° as in Ref. [24]. We 
find a — (89.8^5 4)°. These constraints are in agreement 
with the prediction for a from UTfit [15] and CKMfit- 
ter [14], these predictions are (92.9±5.7)° and (100^7 j)°, 
respectively. 




FIG. 10; Confidence level on a obtained from the flavor SU (3) 
approach. The horizontal dashed lines correspond to the 68% 
(top) and 90% (bottom) CL intervals. The constraint shown 
in this figure uses some of the same inputs as the SU (2) isospin 
constraint of Fig. 9, and the two can not be averaged. 
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XI. CONCLUSIONS 

We report the measurement of the branching fraction, 
fi, and CP violation parameters, Siong and Ciong, for 
the decay P'^P~ using a data sample of (383.6 ± 

4.2) X 10^ SB pairs. We obtain the following results: 

B{B°^p+p-) = (25.5±2.1(stat)+3;g(syst))xl0-^ 

/l = 0.992 ±0.024(stat)+°;°^3(syst), 

^long = -0.17±0.20(stat);°;°^(syst), 
Ciong = 0.01±0.15(stat)±0.06(syst). 

Using these results, and experimental knowledge of the 
other B ^ pp final states, wc perform an isospin analysis 
to obtain a measurement of the CKM angle a. The con- 
fidence level distribution for a is shown in Fig. 9. The so- 
lution obtained that is compatible with the results of SM- 
based fits of existing data is a = [73.1, 117.0]° at 68% CL. 
The uncertainty on the measurement of a from the 
isospin analysis is dominated by penguin pollution. If 
one uses the flavor SU (3) approach described in the 
text to constrain a, one obtains the constraint a = 
[83.3, 105.8]° at 68% CL. 
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APPENDIX: INTERFERENCE OF B° p+ p' 
WITH OTHER B° A-k MODES 

An extensive study of the B backgrounds, and the as- 
sociated systematic error, has been presented earlier in 



this paper. A number of B backgrounds decay into the 
same final state as the signal. These are: 

• B° ^ afTT^.a^ -> p±7rO, 

• ^ p^TT^TT^, non resonant, 

• ^ 7r"'"7r~7r°7r°, non resonant. 

The systematic error associated with these modes was 
evaluated by propagating the uncertainty on the branch- 
ing ratio and CP asymmetries to the; final measure- 
ments, with the appropriate description of acceptance, 
SCF, resolutions and other reconstruction biases (see Sec- 
tion IX B). The likelihood described in Section VII does 
not account for possible interference between 5° P~ 
and other B^ — > 47r final states. 

Acceptance, SCF, and other reconstruction biases are 
not taken into account in this study, and instead, the 
measurements are averaged in a perfect region of 
the phase space referred to as the "p"'"p~ band" , on which 
the analysis technique is assumed to have a uniform sen- 
sitivity. This "p^p~ band" corresponds to the kine- 
matic selection of p^ (see Section IIIB): 0.5 < m7r±7rO < 
l.OGeV/c^ and -0.9 < cos 6»i < 0.98. 

The remainder of this Appendix describes the final 
state wave function, decay amplitudes, effective CP 
asymmetries and the estimate of the systematic uncer- 
tainty from neglecting interference between the signal 
and other B^ Aw final states. 



1. The final state wave functions 

This section summarizes the kinematic dependence of 
the decay amplitudes. 

a. Lorentz invariant phase space 

The four-particle final state can be described com- 
pletely in terms of the following five variables: 

• mi and 1712: the masses of the tt+tt*^ and 7r~7r'^ pairs 

• cos^i and cos ^2: the angle between the tt^tt^ pair 
line of flight (as seen in the B frame) and the in 
the pion pair frame. 

• (p: the angle between the two planes defined by the 
7r+7r° and 'jT~n^ pairs in the B frame. 

We write each phase space dependent quantity rf$4 as a 
function of these five variables. These are: 

oc ^^^^!^dmldmldcos9idcos02d(t), 

Mb TOi 1712 

(34) 

where ^ denotes the momentum of the pion pair in 

the B frame and Mb is the mass of the B (and similarly 
for the other quantities in this expression). 
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b. Kinematic dependence of p 



d. Kinematic dependence of the non-resonant contributions 



We assume that the longitudinal polarization of p~^p~ 
is dominant, so the kinematic dependence of the p~^p~ 
wave function is: 



The kinematic dependence of the non-resonant contri- 
butions can be written as: 



/(p+p-) oc BWimi)BW{m2)Y^''{9i)Yl>{e2) 
oc BW{mi)BW{m2) cos(6'i) cos(6'2), 



(35) 



where BW denotes a Breit-Wigner function (the simple 
non-relativistic form is assumed), and Ij™ are the spher- 
ical harmonics. 



c. Kinematic dependence of am 



In order to write the wave function of tt 



i^TT^, wc need to 
define the following three variables, which are functions 
of the variables found in part la of this Appendix: 

• Op-, the angle between the 7r°and the p line of flight 
(as seen from the oi) in the p frame, 

• 6ai- the angle between the p and Oi line of flight 

(as seen from the B) in the ai frame, 

• (j)ai : the angle between the two decay planes of the 
oi and p mesons, in the ai frame. 

A schematic view of B° ^ 
given in Fig. 11. 



afw^ showing Op and 6ai is 




FIG. 11: Definition of the variables 9p and 9ai entering the 
oiTT kinematic dependence. 

The ai has the possibility to decay both into an S wave 
and a D wave. We assume that the S wave contribution 
is dominant, so the kinematic dependence of ain is: 

fiaiw) (xBW{ai)BW{p) 



(36) 



where the rf^ ^, are rotation matrices. A non-relativistic 

Breit-Wigner is assumed for BW{p), and a relativistic 
form is used for BW{ai). 



/(pT^±^°) (X BW{pi)Y,°{9i), 
for p^7r^7r°, and: 

/(7r'''7r~7r°7r°) = constant, 



(37) 



(38) 



for the non-resonant 4-pion final state. 

We note that the spherical harmonics entering the 
p+p" amplitudes found in Eq. (35) differ from those 
found in these two amplitudes: this means that these 
wave functions arc orthogonal, and after the integra- 
tion over the helicity angles, the interfering terms 
/(p+P")/*(p'^7r±7r°) and /(p+p-)/*(7r+7r-7r°7r°) will 
cancel completely if the integration is done over the whole 
phase space. The "p+ p^ band" used in this analysis is 
slightly asymmetric in helicity angle, and as a result we 
can expect a small deviation from zero for these integrals. 



e. Symmetrization of the wave functions 

The final state is made of four bosons, two of which 
arc identical (the two tt"). Thus each of the above wave 
functions needs to be symmetrized, that is to say that 
the amplitude: 



(39) 



where p denotes the four-momentum of one of the pions, 
needs to be replaced by: 

f{P-K+,Plo,P„-,plo) = -^i^f{p„+,plo,p„-,plo) 

+/(P^+,P^o,p^-,pio)|.(40) 

In the rest of this study, each of the wave functions are 
symmetrized by meson exchange. 

2. Decay amplitudes 

For a particular contribution B° x, the total decay 
amplitude A{B^ x) is the product of a phase space 
independent amplitude A{B'^ x) and one of the above 
kinematic dependences f{x): 

^(B° ^x) = A(B0 ^ x)f{x). (41) 

This section describes how we calculate A{B^ x). 



22 



a. Amplitude of B ^ p 

For B — > p^p~ , two amplitudes need to be known: 

A{B^ p'^ p") and A{B^ p^p~), corresponding to 4 
degrees of freedom. We ctioose the following parameter- 
ization: 

• The branching ratio B{B — > p'^ p~), 

• The CP- violating asymmetries Ciong(p^P~) and 

S'long(P"'"P"), 

• An overall phase 5{p'^ p~) which may be set to zero 
in the appropriate phase convention. 

The parameters B{B p'^p~), Ciong{p'^p~) and 

Siongip~^p~) are set to the values reported here (see Sec- 
tion XI). 



6. Amplitude of B ^ oitt 

Four amplitudes describe the B ^ aiw contribu- 
tions, corresponding to a total of 8 degrees of freedom: 
A{B° a+TT"), A{B° a^TT+), A{B° a^Tr") and 
A{B^ — > aj~7r"'"). We choose a parameterization similar 

to that of the time dependent aivr analysis [67], for which 
parameters have direct physical meaning: 

• B{B'^ oitt) is the total branching ratio. 



C(ai7r) + AC(ai7r) = C+{anr) 
where Acp(a+7r-) = f 



l+\Xcp(a+n-)\^ 



_ l-|A(7p(oi Tr"^ 



Cianr) - AC(ai7r) = C-{anr) 

where Xcp{a^TT+) = ^ ^[laZl'-ll] ^ C(ai7r) is a di- 
rect CP-violating asymmetry and AC{anr) is CP- 
conserving, 



<S'(ai7r) -|- AS{aiTr) = S'+(ai7r) 



2Tm{Acp(a^7r~)} 
l+|Acj:.(a+7r-)|2 ' 



S{anT) is a CP-violating asymmetry in the inter- 
ference between decay and mixing, and AS{aiTr) is 
CP conserving. 



Acp{aiTr) 



\A+-f + \A+-\^-\A- + \^-\A-+\'' 



\A+-\^ + \A+~\^-\A- + \^-\A~^\ 



where 



conjugate. These are also direct CP violating asym- 
metries, 

• Two overall phases of the final states S{af'jT~) and 

For aiTT we assume that all of these parameters are 
constant over the whole of phase space. 

3. Effective CP asymmetries 

Considering the two final states p'^p" and ain, the 
phase-space dependent four-pion decay amplitude is: 

AiB° 47r) = AiB° -> p+p-) + AiB^' a+Tr") 
+ A{B° a^Tr+), 
= A{B^ ^p+p-)f{p+p-) 

+ A{B° ^ at7r-)f{at^-) 
+ A{B° ^a^7r+)f{a^7r+), 

with a similar expression for A{B^ Aw) . Similar equa- 
tions exist for pTrw and 47r final states. 

We can define the phase-space dependent quantity: 



A47r(^4) = 



q A{B° 47r) 
p A{BO An) ' 



(42) 



where A4,r(^4) is a function of mi, 7712, cos^i, cos02, and 
Thus define the phase-space dependent CP-violating 
asymmetries: 



C4x($4) = 

54.($4) = 



1- |A4.($4)P 
1 + |A4.($4)P' 

2Tm(A4^(^4)) 
1 + \X4A^aW'' 



(43) 

(44) 



A^^ = A{B° afw^), and A is its complex 



that are analogous to Eq. (8) and Eq. (9). 

In principle it is possible to measure the value of 
C47r($4) and SiTri^i) at each point of phase space, how- 
ever this would be experimentally challenging. Instead, 
this measurement uses integrated values Cgg and Ses, 
which can be seen as the averages of C and S across the 
previously defined ''p'^p~ band" , weighted over the num- 
ber of events found at each location of the phase space 
which corresponds to: 
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^ /pp band g4.(^4)(|.4(5° ^ 47r)|^ + \A{B^ ^ 47r)p)d^4 

/pp ba„d(l-4(^° - 4^)P + \AB° - 47r)P)d$4 ' 
^ ^ /pp band ^4.($4)(|^(i3" ^ 47r) ^ + lAjB^ A^r)d^, 

/pp band(M(^° - 47r)P + |^(i?0 ^ 47r)|2)d$4 



We can also define the values of C and S that would be 

observed in the absence of interference effects, denoted 
by Cjr(^4) and Sf^{^4), and similarly the integrated 
values across the "p+p- band" C^^ and S^^^. Finally, 
for the evaluation of the systematic error, we are really 
interested in the biases on C and S purely due to inter- 
ference which are the two differences 6C = Ces — 
and 6S = SeS — Sfg. We adopt the same strategy for the 
branching fraction B{p'^p~) by defining 6B = BeS — B^g 
in a similar way. 



elude that interference between signal and other Aw final 
states has a negligible impact on the measurement. 



TABLE VI: Systematic error coming from interference effects 
due to three contaminating amplitudes. 



Contribution 


CT{C) CT{S) ^ 


aiTT 

pTTTT non resonant 

4n non resonant 


0.009 0.009 2.5% 
0.010 0.010 1.1% 

0.002 0.002 0.3% 


Total 


0.014 0.014 2.7% 



4. The systematic error 



As the 8 parameters describing the amplitude 
A{B°{B°) am) are not yet measured, to estimate the 
systematic error on p'^p~ parameters, we need to propa- 
gate the uncertainty on these 8 parameters. For this, we 
use a MC technique. For each experiment, we generate 
the CP parameters C(ai7r), AC(ai7r), ^(aiTr), AS'(ai7r), 
Acp{a.iT^), the phases 5{a'^'K^) and 5(a^7r+), and the 
branching fraction B{aiT:). The values used to generate 
the aiTT contribution have Gaussian PDFs for measured 
quantities (with mean and width corresponding to the 
measurement), and uniform PDFs for unmeasured quan- 
tities (sampled over the allowed values). Then, we esti- 
mate the value oi 6C and 6S for each generated value of 
the aiTT parameters. The systematic error are taken as 
the RMS of the SC, SS and SB distributions. 

For this study, the quantity B{ain) is fixed to the value 
measured in [51], ^(5° ^ of 7r-,a5^ ^ p+7r°) = (16.6 ± 
2.4) X 10~^ and the CP parameters are obtained from [67]: 
dam) = -0.10±0.17, AC(ai7r) = 0.26±0.17, S{anT) = 
0.37 ± 22, A5(ai7r) = -0.14 ± 0.22 and Acpiam) = 
-0.07 ±0.07. 

We estimate the effect of interference between signal 
and either p^tt^tt^ or 7r"'"7r~7r''7r'' in a similar way, where 
we assume the branching fraction for these two non- 
resonant backgrounds is (5.7 ± 13.0) x 10~^. This value 
is obtained by modeling the p^n^Tv'^ component in the 
likelihood fit and by fitting a free yield to the data, and 
we assume that p^Tr^Tr*^ and 7r+7r~7r'^7r'^ have the same 
branching fraction. 

We summarize the systematic uncertainty obtained for 
the different 47r final states considered here in Table VI. 
The total uncertainty on B is used to calculate the quoted 
systematic uncertainty on from this source. We con- 
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